Compressing nearly hard sphere fluids increases glass fragility 
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Abstract. - We use molecular dynamics to investigate the glass transition occuring at large 
volume fraction, ip, and low temperature, T, in assemblies of soft repulsive particles. We find that 
equilibrium dynamics in the [ip, T) plane obey a form of dynamic scaling in the proximity of a 
critical point at T = and Lp — ifo, which should correspond to the ideal glass transition of hard 
spheres. This glass point, 'point G', is distinct from athermal jamming thresholds. A remarkable 
consequence of scaling behaviour is that the dynamics at fixed ip passes smoothly from that of a 
strong glass to that of a very fragile glass as increases beyond po. Correlations between fragility 
and various physical properties are explored. 



Introduction. — Structureless, hard, frictionless par- 
ticles pass from a mobile to an immobile state with in- 
, creasing density [1] . Interacting particles and molecules 
' in a glass-forming material also pass from a mobile fluid 
, state to an immobile glassy state as temperature is re- 
' duced [2] . Much research in the last decade has been 
, devoted to extracting a common geometric essence from 

■ these two classes of phenomena. One line of research ap- 
, proaches the threshold of immobilization or jamming via 
' processes unrelated to thermal equilibrium [3-6] . Connec- 
, tions between these jamming transitions and those seen 
' at positive temperature are suggested [7] , but remain un- 
, clear. In particular, the threshold density for jamming 

■ and that for immobilization at non-zero temperature are 
, considered identical by some researchers [8,9], distinct for 

others [10-12], ill-defined by some others [3,13]. Direct 
measurements are not conclusive, because the location 
of the glass transition relies upon fitting and extrapola- 
tion [8, 10], while the jamming transition is not uniquely 
defined [3,4]. 

The notion that temperature and density should have 
analogous effects on the glass transition has a long his- 
tory [14,15], although quantitative evidence supporting 
these analogies is limited. Recent work studying the effect 
of pressure on the glass transition showed that the dynam- 
ics of glass-formers is little affected by increasing the den- 
sity, since a simple rescaling procedure collapses a broad 



range of dynamic data [16, 17]. This finding directly im- 
plies that the (isochoric) fragility [18] of most glass-formers 
is independent of density, at least in the range currently 
explored by experiments. A second conclusion is that the 
glass transition of molecular systems is mostly controlled 
by temperature, suggesting that the density-driven glass 
transition of hard spheres might have a different nature. 

In this article we study the relative influence of density 
and temperature on the glass transition using a model 
of soft repulsive particles [19]. In the zero-temperature 
limit, the model is equivalent to density-controlled hard 
spheres, while it resembles thermally driven dense fluids 
at large density and finite temperature. Another moti- 
vation to use compressible particles is to access densities 
beyond the hard sphere critical density for kinetic arrest 
that are unreachable with the hard sphere potential. This 
approach has proven useful in the context of athermal jam- 
ming [4-6], but was not extended to thermal equilibrium 
before. At equilibrium, issues related to the possible pro- 
tocol dependence of the results [20] do not arise. 

Using computer simulations we have studied the equi- 
librium dynamics of a three-dimensional assembly of soft 
repulsive particles varying the volume fraction, tp, and 
temperature, T (see fig. 1). We have discovered a sim- 
ple connection between density and temperature effects, 
from which several interesting results are deduced. Equi- 
librium dynamics obey critical scaling in the proximity of 
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Fig. 1: Volume fraction, ip, and temperature, T, phase dia- 
gram of clastic spheres with investigated state points. Four 
iso-relaxation time (dashed) lines are displayed, showing that 
Ta increases by increasing ip or decreasing T. Scaling permits 
accurate determination of the glass point G at T = and 
ipo = 0.635. The glass transition (full) line is obtained assum- 
ing a specific form of the scaling function in eq. (4). 



a glass critical point, which wc call 'point G", at T = and 
If = ipo, see fig. 1. By approaching point G along various 
equilibrium routes, we can determine the functional form 
of the dynamics and the location of the singularity with 
much less uncertainty than is usually possible in glasses. 
In particular, we confirm the non-trivial density depen- 
dence of the relaxation time suggested by recent experi- 
ments on colloidal hard spheres [10]. Scaling also implies 
that the evolution of the dynamics with temperature at 
fixed (f passes smoothly from that of a strong glass to that 
of a very fragile glass as (p increases beyond (po. Compared 
to previous numerical glass models with tunable fragility, 
the present model does not require changing the compo- 
sition of the liquid [21], or the curvature of space [22], 
fragility varies over a much broader range [23], and pro- 
vides a new conceptual way to tune fragility. 

Model and methods. — We use Molecular Dynamics 

simulations [24] to study a system composed of N particles 
enclosed in a periodic cube of linear size L and interact- 
ing through a pair- wise potential: V{rij) = e(l — Vij/aij)"^ 



for rij < aij, V{rij) = otherwise. The interparticle 
distance is = |rj — rj| and aij = {ai + aj)/2, where 
and fTj are the position and diameter of particle i, re- 
spectively. Wc use system sizes between 500 and 8000 
particles, and report results for TV = 1000, for which no 
finite size effects are detected, within numerical accuracy. 
We prevent crystallisation by using a 50:50 binary mixture 
of spheres of diameter ratio 1.4 [4]. The volume fraction 
is if = j^i'^ + 1.43) 0.98/9, with L expressed in units 
of the small particle diameter and p = N/L the number 
density. Up to volume fraction ip = 0.846 we detect no 
sign of crystallization at all studied temperatures; above 
p = 0.924 there was evidence of incipient crystallization 
at the lowest temperatures. However, these crystalization 



effects occur well away from the region of interest around 
Po = 0.635. Wc use e as the energy unit, and \J a\l t as 
time unit, masses are set to unity. All dynamical results 
are obtained at thermal equilibrium, which has been care- 
fully controlled. When temperature is low and density is 
large, we are not able to thermalize. Crystallization and 
equilibrium issues determine the boundaries of the region 
investigated in the phase diagram of fig. 1. 

Dynamic scaling at thermal equilibrium. In 
fig. 2-a we report the evolution of the averaged relaxation 
timescale, ra((/5, T) for all investigated state points. We 
quantify the microscopic dynamics through the self-part 
of the intermediate scattering function: 
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where q is the scattering vector (g ~ 6.1, close to the 
first diffraction peak) and brackets indicate a thermal av- 
erage. We define Tq by -F,((j, To) = e~^, and we arbitrarily 
choose to restrict the average in (1) to the Ni, big parti- 
cles. Data for Tq, are normalized by 1/Vt, which is equiv- 
alent to renormalizing times by the 'thermal' time 

O f MB ' 

where vmb is the first moment of the Maxwell-Boltzmann 
distribution. Therefore, in the T ^ limit where par- 
ticle overlaps are energetically disfavoured, the dynamics 
of the elastic spheres coincide with that of hard spheres 
thermalized at T = 1. We have verified this equivalence 
quantitatively by comparing our numerical results to the 
hard sphere studies presented in Ref. [10]. 

The temperature evolution of Tq, exhibits two qualita- 
tively distinct regimes. For dilute systems, Ta first in- 
creases when T decreases, but saturates when T ^ at 
a finite value corresponding to T^^{<p), the relaxation time 
of the hard sphere fluid. For dense systems, Ta seems to 
increase with no saturation, and we cannot reach equilib- 
rium when T becomes too small and the system becomes 
an unequilibrated glass (see fig. 1). The frontier between 
these two regimes must correspond to p = po, the density 
at which T^^{p) diverges. By analogy with the Jamming 
transition occuring at point J along the p>-axis for ather- 
mal systems of spheres [4], we call 'point G" the location 
of the Glass transition at po along the p-axis for thermal 
systems. 

We now formulate the hypothesis that dynamics in the 
proximity of point G obey scaling behaviour. We surmise 
that elastic spheres at T > behave, in a statistical sense, 
similarly to hard spheres with a reduced diameter [14]. 
Physically, we assume that at low T, overlaps between 
particles are very small, and so are corrections to hard 
sphere behaviour. Below, we justify why an accurate re- 
lation between density and temperature is: 



(2) 



where a > is a numerical prefactor, /i > a charac- 
teristic exponent, and pes{p,T) < p an effective volume 
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Fig. 2: (a) Relaxation timcscales for all investigated state 
points. We rescale Tq, by 1/Vt so that the T ^ hmit coin- 
cides with hard spheres thcrmalized at T = 1. Not all volume 
fractions are labelled, (b) Collapse of the ifi < <po and tp > ipo 
data along the bottom and top branches, respectively using 
eq. (4) and imposing /i = 1.3 from considering the potential 
energy, (c) Arrhcnius plot of the tp > ifo data, using the defi- 
nition Toci'fi,Tg) = 5 ■ 10'^. 



fraction for the elastic sphere system at volume fraction ip 
and temperature T. Furthermore, we build upon a recent 
analysis of the dynamics of colloidal hard spheres [10] and 
assume an exponential divergence for t^^: 



T^^if) ~ exp 



A 



(3) 



where (5 « 2 and ipo ~ 0.637 [10]. Although the exponen- 



tial divergence of t^^ is unambiguous from hard sphere 
studies [10], the values of 5 and po remain subject to large 
uncertainty because the divergence must be extrapolated 
along a single path (increasing (/? at T = 0), relatively far 
away from ipQ. 

We analyze the results of fig. 2-a and approach point G 
from multiple paths in the {<p,T) plane to establish the 
robustness of eq. (3). Combining (2)-(3) we suggest: 



T'a{'P,T) ~ exp 



A 



(4) 



where the scaling functions F±(x) apply to densities 

above/below ipo. We expect therefore that F^{x oo) 
1 to recover the hard sphere fiuid limit, eq. (3), when 
T — > and < ifo- Similarly, F+{x — > oo) — > oo. 
Moreover, continuity of Tq, at finite T and ip = (po implies 
F_(x ^ 0) ~ F+(x ^ 0) ~ x"^^/^. Dynamic scaling was 
recently observed for athermal jamming transitions [5,6], 
but the nature of the critical density and hard sphere di- 
vergence (algebraic instead of exponential) were different 
from eq. (4), while no physical interpretation of scaling in 
terms of an effective hard sphere behaviour was offered. 

The proposed scaling behaviour is confirmed in fig. 2-b 
for data in the range <p> G [0.567,0.736]. To obtain the 
scaling plot, we fix = 1.3 (see below), and use ipo and 
6 as free parameters to collapse |(po — 'pI'^logTa against 
\ipo — 'p>\'^^'^/T. The best collapse is shown, but good results 
are obtained for nearby values of ipo and 5, yielding error 
bars: 

tpo = 0.635 ± 0.005, ^ = 2.2 ± 0.2. (5) 

Outside this range, the collapse quickly deteriorates. Note 

that 6=1, often used to describe hard sphere data [8], 
is inconsistent with our results. Our scaling analysis thus 
lends crucial support to the conclusions of [10]. Of course, 
we cannot exclude that a different dynamic regime can be 
entered when relaxation timescales beyond reach of our 
numerical capabilities are added to the analysis. 

Glass fragility. — The scaling in eq. (4) predicts 

the temperature dependence of Tq, at ipQi Ta{ipQ,T) ~ 
exp(A/r'"'/2). Since iiS/2 « 1.43, this divergence is 
slightly stronger than, but not very different from, the 
simple Arrhcnius behaviour observed for 'strong' glass- 
formers [18]. The divergence of the scaling function F+{x) 
for large argument moreover implies that the temperature 
dependence of Tq for (p > ipa becomes steeper, making 
the materials increasingly 'fragile' [18]. This is vividly 
demonstrated in fig. 2-c, where we conventionally rescale 
T by Tg defined as the temperature where Tq, reaches an 
arbitrary value [18], logiQ Ta{ip,Tg) = Xg. Such a large 
change of fragility was not reported in a particle model 
before [21 23]. Here, it directly results from the interplay 
between (p and T. We quantify fragility by the steepness 
index [18]: 



m = 



9 log 



10 'a 



d{Tg/T) 



(6) 
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Fig. 3: (a) dependence of the steepness index, eq. (6), mea- 
sured from the data in fig. 2-c. The line is the hnear fit pre- 
dicted in eq. (7). (b) The stretching exponent describing the 
time decay of Fs {q, t) for both species coincides with the one 
for hard spheres in the scaling regime of fig. 2-b, uncorrelated 
to fragihty. (c) Temperature dependence of the specific heat at 
ifi = 0.736 obtained by heating the glass slowly. The full line fits 
the T-dependence of cv (T) below Tg , allowing the definition of 
Acv, as shown by the vertical arrow, (d) The (^-dependence of 
Acv (symbols) is well described by a linear increase (full line), 
closely tracking the fragility in (a). 

which increases steadily when (p increases beyond (po, see 
fi-g. 3-a. The linear variation can be rationalized realizing 
that an approximate expression for in eq. (4) can 

be obtained by pushing further our correspondance be- 
tween elastic and hard spheres, and assuming: Ta{f, T) « 
T^[ipefiiip,T)]. This gives F+{x) w (aa;-''/^ - 1)"^ from 
which we get: 

m{ipo + ^if) ^ mo{l + aAip), (7) 

where mo = iJ,SXg/2 and a = a{Xg/A)^^^ . Figure 3-a 
shows that this predicted linear behaviour is accurately 
obeyed over a seven- fold range of fragihty to. A broader 
range of to would be obtained for Xg corresponding to the 
laboratory Tg. Multiplying for instance Xg by a factor 
3 (from Xg = 5 to 15) in Eq. (7) yields to G [18,150], 
quite close to the experimental spectrum for the fragility 
of glass- formers [18]. 

Additionally, the explicit, but approximate, expression 
for F+{x) predicts the location of a glass line, To{ip) = 
[{ip — ipo)/a^^^]'^^^, shown in fig. 1, in the spirit of [15]. 
Unfortunately, our numerical data alone do not allow us to 
determine whether {x) indeed diverges for a finite value 
of its argument. Therefore, while the location of point G 
is very much constrained by our data, the existence of a 
finite temperature singularity for ip > (po remains open. 



Different glasses possess different fragilities, but this 
variability is not satisfactorily understood. Instead, ex- 
perimentalists often correlate the properties of a glass to 
its fragility [18]. Having a model with tunable fragility 
allows theoretical exploration of the part of these corre- 
lations that is not due to the variety of structures (e.g. 
networks, polymer chains) observed in real glass-formers. 
We have studied two popular correlations: the stretching 
of time correlation functions and the specific heat jump at 
Tg. We fit the alpha-relaxation in Fs{q, t) for both species 
with a stretched exponential form: Fg ~ e~*^*/^^ . Wc find 
that /3 is very weakly T-dependent, and report its value 
for different p in fig. 3-b. For (p within the scaling regime 
of fig. 2-b, the observed /3's for small and big particles are 
independent of compression, and hence of fragility, but 
are consistent with those of hard spheres observed pre- 
viously [10]. While some degree of correlation is often 
reported in experiments [2], it is not very strong when 
data for polymers are discarded from the analysis [25] , in 
agreement with our findings. 

Wc measure the potential energy, Vpot(<^,T') = 
(X)i<j specific heat, cy = dVpot/dT. 

Mimicking experiments, we obtain the jump in cy at the 
glass transition by slowly cooling the system, at constant 
ip, down to T = before re-heating at the same rate. For 
all ip, we use very slow rates, rj 3.10"''. The typi- 

cal behaviour of cy(T) upon heating is shown in fig. 3-c. 
As in experiments, a peak is observed when equilibrium 
is recovered. We estimate Acy as the difference between 
the peak height and the glass specific heat, although dif- 
ferent estimates yield qualitatively similar results. The 
(/3-dependence of Acy is shown in fig. 3-d, together with a 
linear fit. As suggested by experiments [2, 18], we find a 
linear relation between fragility and specific heat. 

To qualitatively explain this observation we discuss the 
behaviour of Vpot ■ For a very dilute system, Vpot decreases 
rapidly as T decreases. A collision-based analysis shows 
that V^ot ~ T^/^: the energy decreases faster than linearly 
and vanishes at T = 0. Correspondingly, the 'excess' en- 
ergy above the hard sphere ground state (V^ot = 0) is 
small, explaining the smallness of Acy at small volume 
fraction. For very dense systems, the energy decreases 
linearly at low T, and increases with (p since particles over- 
lap more upon compression. Therefore the larger ip, the 
larger cy in the fluid phase, and the larger is Acy, as ob- 
served in fig. 3-d. These considerations suggest that fur- 
ther analysis of the potential enery landscape properties 
of the present system and comparison with model land- 
scapes [26,27] could be very valuable. 

In a broad density range encompassing (po, simulations 
suggest Vpot ~ T^, with /x w 1.3, a value intermediate 
between the dilute and dense limits. Due to the harmonic 
nature of the potential, y/Vp^ represents the average over- 
lap between interacting particles, and suggests a way to 
quantitatively estimate ipcs in eq. (2), and to justify the 
form of the scaling variable in eq. (4) . These energetic con- 
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siderations are a suggestive physical interpretation, rather 
than a rigorous derivation, of the success of eq. (2) at col- 
lapsing data in fig. 2-c. 

Non-equilibrium jamming at T = 0. Finally, 
the low-T, non-equilibrium behaviour of the energy den- 
sity during slow annealing is also informative. For ^p < 

Lp* 0.662, V^ox_(lp,T ^ 0) = 0; while Vpot remains finite 
above Therefore, the nature of the T = glasses pro- 
duced by slow annealing at volume fractions above point 
G changes at ip* where a 'jamming' transition similar to 
the one described in [4] occurs. Note that kp* is larger 
than ipj f« 0.648, the critical density for jamming deter- 
mined in [4] for the same system, because our glasses have 
been annealed-with no trace of demixing or crystalliza- 
tion [3, 13]. We emphasize that these jamming transitions 
occur at different densities when different protocols are 
used [13], while the equilibrium glass transition at an- 
alyzed in this work is uniquely defined by the behaviour of 
the equilibrium relaxation time, so that both transitions 
should be considered as distinct phenomena [12]. 

Discussion and conclusion. — The dynamic scaling 
behaviour in cq. (4) describing the interplay between den- 
sity and temperature for soft repulsive particles is in stark 
contrast with the results obtained for soft spheres inter- 
acting with a pure inverse power law potential [16, 17], 
because the interaction potential of elastic spheres has 
a finite cutoff (the particle diameter). Therefore, elastic 
spheres reduce to hard spheres in the small-T limit, which 
is crucial to obtain a change in fragility. We believe our 
results would generically carry out for repulsive potentials 
with a finite range. 

A second interesting feature of this study is that the lo- 
cation of a divergence for the equilibrium relaxation time 
is extrapolated from a demanding scaling procedure where 
a large set of independent data is used to locate a critical 
density. To the best of our knowledge, such an analysis us- 
ing two control parameters has no counterpart in the glass 
transition literature. In particular, we believe it yields a 
rather accurate determination of the critical density cpo 
where the equilibrium relaxation time of the hard sphere 
fluid extrapolates to oo, but leaves open the existence of a 
finite temperature singularity above (po- Our determina- 
tion of the existence of a genuine divergence of Tq, for hard 
spheres at is based on the analysis of the metastable 
fluid demixing and crystallization being avoided due to 
polydispersity. Whether this ideal glass transition is itself 
avoided [13] due to crystallization or demixing at much 
larger density is another important issue that our data 
leave open. 

In conclusion, our preliminary studies of repulsive elas- 
tic particles at thermal equilibrium suggest that soft re- 
pulsive particles are a promising new tool to gain deeper 
understanding of glass transition phenomena in colloidal 
and molecular systems. Experimentally, this could be di- 
rectly realized using soft colloidal particles. 
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